arXiv:1502.02986vl [hep-ph] 10 Feb 2015 


Kinematic biases on centrality selection of jet events in 
pPb collisions at the LHC 


Nestor Armesto 

Departamento de Fisica de Particulas and IGFAE, 

Universidade de Santiago de Compostela, E-15706 Santiago de Compostela, Calicia-Spain 

Doga Can Giilhan 

Laboratory for Nuclear Science and Department of Physics, 

Massachusetts Institute of Technology (MIT), Cambridge, MA 02139 USA 

Jose Guilherme Milhano 

CENTRA, Instituto Superior Tecnico, 

Universidade de Lisboa, Av. Rovisco Pais, P-1049-001 Lisboa, Portugal 

and 

Physics Department, Theory Unit, CERN, CH-1211 Ceneve 23, Switzerland 


Abstract 

Centrality selection has been observed to have a large effect on jet observables 
in pPb collisions at the Large Hadron Collider, stronger than that predicted by 
the nuclear modification of parton densities. We study to which extent simple 
considerations of energy-momentum conservation between the hard process and 
the underlying event affect jets observables in such collisions. We develop a 
simplistic approach that considers first the production of jets in a pp collision 
as described by PYTHIA. From each pp collision, the value of the energy of the 
parton from the proton participating in the hard scattering is extracted. Then, 
the underlying event is generated simulating a pPb collision through HIJING, 
but with the energy of the proton decreased according to the value extracted 
in the previous step, and both collisions are superimposed. This model is able 
to capture the bulk of the centrality effect for central to semicentral collisions, 
for the two available sets of data: dijets from the CMS Collaboration and single 
jets from the ATLAS Collaboration. As expected, the model fails for peripheral 
collisions where very few nucleons from Pb participate. 
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1. Introduction 


The classification of events according to some measurement of their activity 
(energy, multiplicity, etc.) in some region of phase space, what is commonly 
referred to as centrality selection, is extensively used in the study of nuclear 
collisions. Naively one expects that to events within a given centrality class 
correspond similar conditions, say of energy density or temperature. The de¬ 
pendence of specific observables on centrality should then yield information on 
their sensitivity to the properties of the system (energy density, temperature, 
etc.). Further, centrality should provide a link, albeit in a model dependent 
way, to theoretically well-defined quantities such as the impact parameter of 
the collision. 

The usefulness of centrality selection relies, however, on how robustly it can 
be defined. Specifically, whether the centrality selection (and its link to some 
theoretical quantity as the impact parameter) depends (i) on the centrality 
criterium; and (ii) on the observable under consideration i.e. whether there is 
a correlation between the characteristics of the events under consideration, say 
the requirement of presence of jets, and the centrality criterium. 

While in non-peripheral collisions between heavy ions the classification of 
centrality classes appears to be robust, in proton-lead collisions at the Large 
Hadron Collider (LHC) the situation is far more problematic]^ Whereas for both 
dijets [3] and single jets [4] the minimum bias results are well reproduced by 
standard pQCD with nuclear modification of parton densities [5], the centrality- 
selected results show a strong dependence on centrality that cannot be accom¬ 
modated by the existing ideas on the impact parameter dependence of nuclear 
parton densities lilTlIHl- There is an ongoing discussion on the definition of 
centrality in such asymmetric systems laiin]. 

A numbers of factors can conceivably be at the origin of the observed effects. 
Fluctuations in the number of partipating nucleons are much larger in proton- 
nucleus than in nucleus-nucleus collisions and certainly confound the connection 
of centrality to collision impact parameter. Further, energy-momentum conser¬ 
vation poses significantly more stringent demands in the case of asymmetric 
collisions. Ultimately the problems can be traced back to our lack of under¬ 
standing of the detailed microscopic dynamics underlying soft particle produc¬ 
tion in hadronic and nuclear collisions, and its coupling to hard subprocesses. 
These issues have been discussed in several recent works [mm min]. 

In this short note we address solely, and in a simplistic framework, the role of 
kinematic constraints on the energy-momentum of the proton. We ask to what 
extent the requirement of a hard subprocess restricts the soft particle production 
which will eventually determine the centrality of the evenlj^ 


^The problem also exists in dAu collisions at the Relativistic Heavy Ion Collider (RHIC), 

see e.g. [I][2]. 

^Take for example the production of a dijet pair at pseudorapidity »7dijet = 2 with leading 
jet pj’ = 120 GeV. The minimum momentum fraction taken from the proton is then Xp ~ 0.35. 
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We focus on CMS dijet [3] and ATLAS single-jet [3] results in pPb colli¬ 
sions at 5.02 TeV/nucleon. We find that our simplistic implementation leads to 
large effects in good agreement with data. Importantly the model fails to de¬ 
scribe data when the implementation is evidently deficient, that is when energy- 
momentum constraints on the Pb nucleus are expected to play a role. 

The manuscript is organised as follows: in section we describe the model 
that implements energy-momentum conservation constraints on the proton. 
Then, in sections and we employ the model to analyse the extent to which 
CMS results on dijets and ATLAS results on single jets, respectively, can be 
understood as resulting from kinematic considerations. Finally, in section]^ we 
present our conclusions. 

2. Description of the model 

The model that we employ uses PYTHIA m for the hard scattering and 
HIJING [inilT^ for the underlying event. Each event (our results are based on 
samples of 10® for dijets and 2 • 10® for single jet spectra) is generated according 
to the following procedure: 

1. generate a pp event in PYTHIA with no underlying event (for a pro¬ 
ton beam of Ap = 4 TeV and a Pb beam of Apb = 1.58 TeV/nucleon, 

= 5.02 TeV) with the required characteristics: jets or dijets within 
the experimental kinematic cuts, and extract the momentum fraction Xp 
of the hard parton from the proton participating in the 2 —>■ 2 hard sub¬ 
process. The nuclear modification of parton densities was not taken into 
account because for this kinematics their effect, as discussed previously 
(see also [S]), is much smaller than the one observed in data; 

2. generate a minimum bias pPb event in HIJING, for a proton beam of 
(1 — xp)Ep and a Pb beam of Apb, i.e. we reduce the proton energy for 
the underlying event such that y/s^N = 2-^/(1 — Xp)EpEp\,; 

3. shift both events (that are generated in their respective center of mass) to 
a common frame (the LHC laboratory frame); 

4. the HIJING event, generated as minimum bias, is reweighted according to 
its impact parameter such that scaling in the number of collisions Ncoii, 
expected for hard events, is fulfilled; 

5. the PYTHIA and HIJING events are superimposed, resulting in a full pPb 
event with the weight given in the previous step. 

Note that the momentum fraction xpb of the parton from the Pb participat¬ 
ing in the 2 —> 2 hard subprocess is unchanged. This is clearly a deficiency of 
the model, which should accordingly be expected to fail for small Ncou (when 
the number of participating nucleons from the Pb is small) and for large xpb 
(for the Pb-going side in pseudorapidity). 
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3. CMS dijet results 


In this section we consider the dijet measurements by CMS [5]. We generate 
hard events in PYTHIA with jets reconstructed within |? 7 | < 3 using the anti- 
fc-r sequential recombination algorithm nans] with a distance parameter of 
0.3, considering events required to have a dijet with a leading jet pT,i > 120 
GeV/c, a subleading jet pt ,2 > 30 GeV/c, and an azimuthal distance between 
them A(j)i 2 > Then, for the full PYTHIA+HIJING generated event, the 

transverse energy in the region 4 < \r]\ < 5 is rescaled, to match that measured 
by CMS, by a factor that accounts for detector resolution effects that are not 
corrected by CMS. It is this rescaled energy that is used to classify events in 
centrality. 

In figure we show the comparison of our results for the average dijet pseu¬ 
dorapidity (?7dijet = {vi + ^ 2 )!'^) as a function of the total transverse energy 
deposition in the pseudorapidity range of 4 < I 77 I < 5, compared with data mea¬ 
sured by CMS [3|. The agreement is very good except for the lowest transverse 
energies, corresponding to peripheral collisions where the model is not expected 
to work as discussed in the previous section. In figure]^ we show the ratio of the 
dijet pseudorapidity distribution for events in a a given centrality class over the 
minimum bias distribution, compared to the CMS data. Again, a good overall 
agreement is found. 
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Figure 1: Average dijet pseudorapidity (r/^ijet = ixii + ^2)/2) as a function of the total 
transverse energy deposition in the pseudorapidity range of 4 < |r;| < 5 for PYTHIA+HIJING 
events matched according to Xp value (black squares) overlaid with data measured by CMS 
[3] (red circles). 


CMS [3] also provided data on the average dijet pseudorapidity for a fixed 
energy in the p-going direction, for 4 < ?; < 5, as a function of the energy 
in the Pb-going direction, for —5 < rj < —4. In figure we show the 
comparison of the CMS data with our PYTHIA-I-HIJING results. Once more, 
the model captures the trend of data except for the lowest activity data where 
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Figure 2: Ratios of dijet pseudorapidity (?7dijet = (^i +^2)/2) distributions with a selection on 
total forward energy deposition to the dijet pseudorapidity distribution without 

any requirement on event activity. The calculation in PYTHIA+HIJING is shown by black 
squares and CMS data points are shown by red circles. 


only one or two nucleons from the Pb nucleus contribute and the model - that 
corrects for energy only on the proton side - is clearly deficient due to the neglect 
of energy-momentum constraints for Pb. 

4. ATLAS jet results 

In this section we consider the single-jet measurements by ATLAS [1]. We 
generate hard events in PYTHIA with jets reconstructed using the anti-ZcT se¬ 
quential recombination algorithm [IHIIIS] with a distance parameter of 0.4, in 
the region jTyjet — ?7 cm| < 3. The centrality criterium is, in this case, the to¬ 
tal transverse energy in Pb-going direction within the pseudorapidity range of 
—4.9 < 77 < —3.2, thus less sensitive to energy constraints on the proton. Let 
us note that we use, for the different centrality classes, the number of collisions 
NcoU provided by ATLAS and not the one extracted from HIJING|^ 


^They are (11.94,9.86,8.38,6.934,4.82,2.29) for HIJING and 

(14.57,12.07,10.37,8.94,6.44,2.98) in the ATLAS paper [4], for the 0 - 10%, 10 - 20%, 
20 — 30%, 30 — 40%, 40 — 60% and 60 — 90% centrality classes respectively. 
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Figure 3: (Left) CMS data [3], (Right) PYTHIA+HIJING with Xp matching. Colored markers 
show average dijet pseudorapidity (??dijet = + ^ 2 )/ 2 ) as a function of transverse energy 

deposition in the Pb-going direction, for bins of the transverse energy deposition in the p- 

going direction, El^. The choice of pseudorapidity intervals for E^^ and E!^ are —5 < rj < —4 
and 4 < < 5 respectively. The black line shows average dijet pseudorapidity as a function 

of Ejp’ when the requirements on transverse energy in proton going direction are removed. 


In figure we show the results of the model for the nuclear modification 
factor Rppb of jets as a function of their transverse momentum, for different 
centrality classes, integrated over the whole pseudorapidity region |?7jet — ??cm| < 
3. The effect of the centrality selection becomes evident. 

In figures and we show a comparison of the results of the model with 
ATLAS data [3] for the nuclear modification factor Rppb of jets as a function of 
their transverse momentum in different pseudorapidity bins, for central, semi¬ 
central and peripheral collisions respectively. A good agreement with data is 
found for central collisions that deteriorates with decreasing centrality, until the 
model fails for peripheral collisions as expected from previous discussions. 

5. Conclusions 

In this work we have analyzed the effect of centrality selection on events 
characterised by a hard observable in pPb collisions at the LHC. We have fo¬ 
cused on the influence of the interplay given by simple considerations of energy- 
momentum conservation, between the hard process that is the observable of 
interest - jet production - and the underlying event that provides the central¬ 
ity estimator. We have developed a simplistic model that considers first the 
production of jets in a pp collision as described by PYTHIA. From each pp 
collision, the value of the energy of the parton from the proton participating 
in the hard scattering is extracted. Then, the underlying event is generated 
simulating a pPh collision through HIJING, but with the energy of the proton 
decreased according to the value extracted in the previous step. The full event 
consists in the superposition of the hard and underlying ones. In this way the 


6 








: PVTHlA+HIJING 

: Xp matched 


♦ - 

■ h "ii l<3 i 

_ ' 'jet 'cm' T _ 


<1 

; 0-10% 

i 10-20% 

; 20-30% 


,i 


• 

• 

• 

• 

• 

• 

• 

• 



30-40% 

40-60% 

1 

; 60-90% 

10^ 

P, (Gev/c) 

10^ 

P, (GeV/c) 

10^ 

(GeV/c) 


Figure 4: Rppb calculation in PYTHIA+HIJING with Xp matching are shown for 0 — 10%, 
10 — 20%, 20 — 30%, 30 — 40%, 40 — 60% and 60 — 90% centrality classes. The centrality 
classes are determined according to total transverse energy in Pb-going direction within the 
pseudorapidity range —4.9 < rj < —3.2. 
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Figure 5: RpPb calculation in PYTHIA+HIJING with Xp matching (grey bands) and the 
measured values by ATLAS ([4], red circles) are shown for the 0 — 10% centrality class in 
bins of pseudorapidity in the center-of-mass frame. The centrality classes are determined 
according to total transverse energy in Pb-going direction within the pseudorapidity range 
-4.9 < 77 < -3.2. 


energy constraints on the proton are taken into account, while those on the Pb 
nucleus are not considered. The model is thus expected to fail for peripheral 
collisions where very few nucleons from Pb participate. Note that our model 
is not a dynamical one, see other explanations in [miiaiisiiii], our only aim 
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Figure 6 : RpPb calculation in PYTHIA+HIJING with Xp matching (grey bands) and the 
measured values by ATLAS (HI. blue circles) are shown for the 20 — 30% centrality class 
in bins of pseudorapidity in the center-of-mass frame. The centrality classes are determined 
according to total transverse energy in Pb-going direction within the pseudorapidity range 
-4.9 < 77 < -3.2. 



Figure 7: RpPb calculation in PYTHIA+HIJING with Xp matching (grey bands) and the 
measured values by ATLAS (a. black circles) are shown for the 60 — 90% centrality class 
in bins of pseudorapidity in the center-of-mass frame. The centrality classes are determined 
according to total transverse energy in Pb-going direction within the pseudorapidity range 
-4.9 <ri < -3.2. 


being the study of how important this class of biases may become on existing 
experimental data. 

We have considered two sets of data: dijets from CMS [3] and single jets from 
ATLAS [3] and find that the model is able to capture the bulk of the centrality 
effect for central to semicentral collisions, while it fails - as it should - for pe¬ 
ripheral collisions. We conclude that this simple bias due to energy-momentum 
conservation must be considered before any extraction of other eventual cen¬ 
trality dependencies as those of nuclear parton densities. 

The obvious extension of this work is to include the energy-momentum con- 







































straints on the nucleus side, which we plan to address. In any case, we find that 
data as those discussed in this manuscript, offer most valuable information for 
the development of models for high-energy hadronic and nuclear collisions, as 
they link perturbative computable quantities as jets with those coming from the 
underlying event that can only be addressed through accurate modelling soft 
particle production. The present discussions on centrality definition in asym¬ 
metric collisions make it clear that such models are critically needed. 
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